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Abstract 

We derive concentration inequalities for functions of the empirical 
measure of large random matrices with infinitely divisible entries and, 
in particular, stable ones. We also give concentration results for some 
other functionals of these random matrices, such as the largest eigen- 
value or the largest singular value. 
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1 Introduction and Statements of Results: 

Large random matrices have recently attracted a lot of attention in fields such 
as statistics, mathematical physics or combinatorics (e.g., see Mehta [22], Bai 
and Silverstein [3], Johnstone |16| . Anderson, Guionnet and Zeitouni pQ). For 
various classes of matrix ensembles, the asymptotic behavior of the, properly 
centered and normalized, spectral measure or of the largest eigenvalue is 
understood. Many of these results hold true for matrices with independent 
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entries satisfying some moment conditions (Wigner [32], Tracy and Widom 



[5U] , Soshnikov [SB], Girko [7], Pastur [23], Bai [2J, Gotze and Tikhomirov 



There is relatively little work outside the independent or finite second mo- 
ment assumptions. Let us mention Soshnikov [28] who, using the method of 
determinants, studied the distribution of the largest eigenvalue of Wigner ma- 
trices with entries having heavy tails. (Recall that a real (or complex) Wigner 
matrix is a symmetric (or Hermitian) matrix whose entries M^, 1 < i < N, 
and Mj j, 1 < i < j < N, form two independent families of iid (complex val- 
ued in the Hermitian case) random variables.) In particular, (see [28]), for 
a properly normalized Wigner matrix with entries belonging to the domain 
of attraction of an a-stable law, liniTv^oo ^ N {^max < x) — exp (— x~ a ) (here 
Amax is the largest eigenvalue of such a normalized matrix). Soshnikov and 
Fyodorov [29] further derived results for the largest singular value of K x N 
rectangular random matrices with independent Cauchy entries, showing that 
the largest singular value of such a matrix is of order K 2 N 2 . 

On another front, Guionnet and Zeitouni [9J, gave concentration results 
for functionals of the empirical spectral measure for random matrices whose 
entries are independent and either satisfy a Logarithmic Sobolev inequality 
or are compactly supported. They obtained in that context, the subgaussian 
decay of the tails of the empirical spectral measure when deviating from its 
mean (see also Ledoux [IE])- Our purpose in the present work is to deal 
with matrices whose entries form a general infinitely divisible vector, and in 
particular a stable one. We obtain concentration results for functionals of the 
corresponding empirical spectral measure, allowing for any type of light or 
heavy tails. The methodologies developed here apply as well for the largest 
eigenvalue or for the spectral radius of such random matrices. 

Following the lead of Guionnet and Zeitouni [SJ, let us start by setting 
our notation and framework. 

Let A^ATxAf(C) be the set of N x N Hermitian matrices with complex 
entries, throughout, equipped with the Hilbert-Schmidt norm 



Let / be a real valued function on R. The function / can be viewed as map- 
ping M NxN (C) to M NxN (C). Indeed, for M = (My)i< y < s G M NxN (C), 



N 
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so that M = UDU*, where D is a diagonal matrix, with real entries Ai, Xn, 
and U is a unitary matrix, set 



f(M) = U/(D)U*, /(D) 



/ /(Ai) ■•■ \ 

/(A 2 ) ••• 



••■ f(X N ) J 
Let £r(M) = J2?=i ^i,; be the trace operator on A4nxn(C) and set also 

1 N 

tr N (M) = -J2 M ^- 
i=i 

For a A x A random Hermitian matrix with eigenvalues Ax, A2, Xn, let 
Fn(x) = jr J2iLi l{A t <x} be the corresponding empirical spectral distribu- 
tion function. As well known, if M is a iV x N Hermitian Wigner ma- 
trix with E[Mi,i] = E[M lj2 ] = 0, E[|M lj2 | 2 ] = 1, and E[Mf )]L ] < 00, 
the spectral measure of W/l/y/N converges to the semicircle law: a(dx) = 
V4^1 {lxl < 2} dx/2n (P). 

We study below the tail behavior of either the spectral measure or the 
linear statistic of /(M) for classes of matrices M. Still following Guionnet 
and Zeitouni, we focus on a general random matrix Xa given as follows: 

Xa = ((XA)ij)i<y<jv, Xa = X A , (Xa)ij = -y=AijUJij, 



with {oJi,j)x<i,j<N = + y/-lu! J )x<i ) j<N, = and where 1 < 
i < j < N is & complex valued random variable with law Pjj = P^ + 



-lPjj, 1 < % < j < N, with P^ = So (by the Hermite property). Moreover, 
the matrix A = (-Aij)i<ij<jv is Hermitian with, in most cases, non-random 
complex valued entries uniformly bounded, say, by a. 

Different choices for the entries of A allow to cover various types of ensem- 
bles. For instance, if Uij, 1 < % < j < N, and u^i, 1 < % < N, are iid N(0, 1) 
random variables, taking A^i = v2 and A it j = 1, for 1 < i < j < N 
gives the GOE (Gaussian Orthogonal Ensemble). If u^j,uj^, 1 < i < j < N, 
and uf 1 ^ 1 < i < N, are iid A^O, 1) random variables, taking A iti = 1 and 
Aij = l/y/2, for 1 < i < j < N gives the GUE (Gaussian Unitary Ensem- 
ble) (see [22]). Moreover, if lu^, u(j, 1 < i < j < N, and u^, 1 < i < N, are 
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two independent families of real valued random variables, taking Aij = for 
\i — j\ large and A+j = 1 otherwise, gives band matrices. Proper choices of 
non-random A^j also make it possible to cover Wishart matrices, as seen in 
the later part of this section. In certain instances, A can also be chosen to be 
random, like in the case of diluted matrices, in which case Aij, 1 < i < j ; < N, 
are iid Bernoulli random variables (see Pj). 

On R^ 2 , let be the joint law of the random vector X = (uf\, ui^, lo\ ■), 
1 < i < j < A^, where it is understood that the indices for u; t " are 1 < i < 
N. Let be the corresponding expectation. Denote by fi 1 ^ the empirical 
spectral measure of the eigenvalues of Xa, and further note that 

^/(X A ) = ltr(/(X A )) = / f{x)fil{dx), 
for any bounded Borel function /. For a Lipschitz function / : W d — > M, set 

Lip ~^v \\x-v\\ ' 

where throughout || • || is the Euclidean norm, and where we write / G Lip(c) 
whenever \\f\\ Lip < c. 

Each element M of M.nxn(^-) has a unique collection of eigenvalues A = 
A(M) = (Ai, • ■ • , Ajy) listed in non increasing order according to multiplicity 
in the simplex 

S N = {Ai > ■ ■ ■ > \ N : \i E R, 1 < i < N}, 

where throughout S N is equipped with the Euclidian norm ||A|| = v/X^Li A?. 
It is a classical result sometimes called Lidskii's theorem ([24]). that the 
map .A/1 7v x TV (C) — > S N which associates to each Hermitian matrix its or- 
dered list of real eigenvalues is 1-Lipschitz ([10j. [17]). For a matrix Xa 
under consideration with eigenvalues A(X A ), it is then clear that the map 
if : (wg, u{J)i<i<3<N ^ A(X A ) is Lipschitz, from (R N * , \\ ■ ||) to (S N , || ■ ||), 
with Lipschitz constant bounded by a^j2/N. Moreover, for any real val- 
ued Lipschitz function F on S N with Lipschitz constant H-FHl^, the map 
F o ip is Lipschitz, from (M^ 2 , || ■ ||) to R, with Lipschitz constant at most 
a\\F\\Li P ^2/N. Appropriate choices of F ([17], P) ensure that the maxi- 
mal eigenvalue XmaxO^x) = Ai(X A ), the spectral radius o(Xa) = max \\A 

Ki<N 
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and irjv(/ (Xa)), where / : R — > R is a Lipschitz function, are themselves 
Lipschitz functions with Lipschitz constants at most a^2/N, a^j2/N and 
V / 2 fl 11/ II Li P /N, respectively. These observations (and our results) are also 
valid for the real symmetric matrices, with proper modification of the Lips- 
chitz constants. 

Next, Recall that X is a <i-dimensional infinitely divisible random vector 
without Gaussian component, X ~ ID((3,0,i>), if its characteristic function 
is given by, 

<p x (t) = Ee^ x) 

= exp ji(t,/3) + (e 1 ^ - 1 - i^u)!^) v(du)} , (1.1) 

where t, (3 G R d and v (the Levy measure) is a positive measure on 
B(R d ), the Borel a-field of R d , without atom at the origin, and such that 
L d (l A ||tt|| 2 )z/((iM) < +oo. The vector X has independent components if 
and only if its Levy measure v is supported on the axes of R d and is thus of 
the form: 

d 

v[dx u . . . , dx d ) $o{dxi) . . . 5 (dx k ^ 1 )P k (dx k )8 (dx k+ i) . . . 5 (dx d ), (1.2) 

k=l 

for some one-dimensional Levy measures v k . Moreover, the v k are the same 
for all k — 1, . . . , d, if and only if X has identically distributed components. 

The following proposition gives an estimate on any median (or the mean, 
if it exists) of a Lipschitz function of an infinitely divisible vector X. It 
is used in most of the results presented in this paper. The first part is a 
consequence of Theorem 1 in |13j . while the proof of the second part can be 
obtained as in [T3] . 

Proposition 1.1 Let X = (ufi, iof^ ^u)i<i<j<N ~ ID(/3,0,v) in R^ 2 . Let 

y2 ( x ) = I\\u\\<x \\ u \\ 2iy ( du ), 9 ( x ) = S\\ u \\>x V ( du )> and for any 7 > 0, let 
p-y = inf {x > : < V 2 (x)/x 2 < 7}. Let f G Lip{l), then for any 7 such 
that u(p 7 ) < 1/4, 

(i) any median m(f (X)) of f(X) satisfies 

\m{f(X))-f{0)\ < Gxirf) :=p 7 (V7 + 3fc 7 (l/4)) +£ 7 , 
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(ii) the mean E, N [f(X)} of f(X), if it exists, satisfies 

\E N [f(X)\ -/(0)| < G 2 ( 7 ) :=p 7 (v^+ fc 7 (l/4)) + £ 7 , 

where k 7 (x), x > 0, is the solution, in y, of the equation 

y - (y + 7) in ^1 + ^ = ln ^, 

and where 

( n2 r r A 1/2 

£ 7 = y](( e fe> & )- / {e k ,y)v(dy)+ \ (e k ,y)u(dy)) 

\ /, , v ^<l|y||<i ■/i<lli/ll<P7 y / 

(1.3) 

wit/i ei, e 2 , . . . , e/v2 being the canonical basis ofW N2 . 

Our first result deals with the spectral measure of a Hermitian matrix 
whose entries on and above the diagonal form an infinitely divisible random 
vector with finite exponential moments. Below, for any b > 0, c > 0, let 

Lip h {c) = {/ : R -> R : ll/H^ <c, ||/||oo < &}, 

while for a fixed compact set /C C R, with diameter |/C| = sup |x — y|, let 

x,y€K 

U V k{c) ■= {/ : R - R : 11/11^ < c, C /C}, 

where supp(f) is the support of /. 

Theorem 1.2 Lei X = (u)(\, ixifj, w(j)i<i<j<N be a random vector with joint 
law F N ~ ID(/3,0,v) such that E^e*"*"] < +oo 7 for some t > 0. Let 
T = sup{t > : E^fe* 11 * 11 ] < +00} and let h' 1 be the inverse of 

h(s) = [ \\u\\ (e" M - l)u(du), < s < T. 
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(i) For any compact set K, C R, 

W sup \tr N (f(X A ))-E N [tr N (f(X A j)]\>5 



VGLip K (l) 

for all5>0 such that 5 2 < %y/2a\K,\h (T~) /N 



S\JC\ f 8via|K| . 

<^d exp J / h- 1 (s)ds}, 1.4) 



P N ( sup \tr N (f(K A ))-E N [tr N (f(K A ))]\ >6 



v /eLi P6 (i) 



C((5, b) I f V2aC(S,b) 



<^^exp<J - / h- 1 (s)ds}, (1.5) 



/or a// 5 > sitcA that 5 2 < V2aC{5, b)h(T~)/N, wh 



ere 



C{6,b) = c(^(G 2 W+h{taj) + 



with (^2(7) as in Proposition \1. 11 C a universal constant, and with t 
the solution, in t, of th(t) — J Q h(s)ds — ln(126/5) = 0. 

Remark 1.3 (i) The order of C(5,b) in part (ii) can be made more spe- 
cific. Indeed, it will be clear from the proof of the theorem (see 112. 39) ), 
that for any < t* < T fixed, 

, . fV2a/ln^ [**h(s)ds 



(ii) As seen from the proof (see \2. 39) ). in the statement of the above the- 
orem, (^2(7) can be replaced by K N [||X||] . Now K N [||X||] is of order 
N since, 

N min E N \\XA] <E N \\\X\\\ < N max E N \X?], (1.6) 

3=1,2,...,JV2 Ll J ' i L " " J j=l,2,...,JV 2 1 3i 
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where the Xj, j = 1, 2, . . . , N 2 are the components of X . Actually, an 
estimate more precise than U.6]) is given by a result of Marcus and 
Rosinski JW$ which asserts that ifE[X] = ; then 

~x <E[||X||] <^x , 
where xq is the solution of the equation: 

YM + m =h (1 . 7) 

x 2 x 

where V 2 (x) is as before, while M(x) = L u » >x \\u\\v{du) , x > 0. 

(Hi) As usual, one can easily pass from the meanE N [trN(f)] to any median 
m(tr N (f)) in either or U.5\) . Indeed, for any < 5 < 2b, if 

sup \tr N (f) - m(tr N (f))\ > 5, 
feLi Pb (i) 

there exist a function f G Lip^l) and a median m(tr/v(,/)) of tr^{f), 
such that either tr N (f) — m(trjv(/)) > 5 or tr N (f) — m{tr^{f)) < —5. 
Without loss of generality assuming the former, otherwise dealing with 
the latter with —f, consider the function g(y) = min (d(y, A), 5) /2, 
y G M. N , where A = {tr^{f) < m{tr^{f)} . Clearly g G Lipb(l), 
E N [tr N (g)] < 5/4, and therefore tr N {g) - E N [tr N (g)] > 5/4, which 
indicates that 

r 

sup \tr N {g)-E N [tr N {g))\ > - 

9ELip b {l) 4 



Hence, 



sup \tr N (f) - m(tr N (f)) \ > 5 
feLi Pb (X) 



<F N ( sup \tr N (g)-E N [tr N (g)]\> 5 -). (1.8) 

\geLip b (l) V 

Next, recall (see [BJ, [I?]) that the Wasserstein distance between any two 
probability measures Hi and /i 2 on M is defined by 



d w {fi 1 ,fi 2 )= sup / /d/ii - / 

f£Lip b (l) Jr Jm 
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(1.9) 



Hence, Theorem II .21 actually gives a concentration result, with respect to the 
Wasserstein distance, for the empirical spectral measure when it deviates 
from its mean E^/t^]- 

As in [9], we can also obtain a concentration result for the distance be- 
tween any particular probability measure and the empirical spectral measure. 

Proposition 1.4 Let X = (^^,^^,^fj)i<i<j<jv be a random vector with 
joint law F N ~ ID(f3, 0, v) such that E N [e* l|x|i ] < +oo ; for some t > 0. Let 
T = sup{t > : E^JV^"] < +00} and let h^ 1 be the inverse of h(s) = 
f RN 2 ||u||(e s """ — l)u(du), < s < T. Then, for any probability measure /i, 

V N (dwifi 1 !^) - ® N [dw(P-A,»)] > s ) < ex P I " h- 1 (s)ds\, (1.10) 

for allO < 5 < V^ah (T~) /N. 

Of particular importance is the case of an infinitely divisible vector having 
boundedly supported Levy measure. We then have: 

Corollary 1.5 Let X = (u;^, ui^, &lj)i<i<j<N be a random vector with joint 
law F N ~ ID ((3, 0, v) such that v has bounded support. Let R = inf{r > : 
v{x : \\x\\ > r) = 0}, let V 2 ( = V 2 (R)) = j RN 2 ||w|| 2 z/(<i-u) ; and let 

£(x) = (l + x) ln(l + x) -x, 

x > 0. 

(i) For any 5 > 0, 



W sup \tr N (f(X A ))-E N [tr N (f(X A ))]\>5 
v /eLi P6 (i) 



with 6*2(7) as i n Proposition li.il C a universal constant, and to the 
solution, in t, of 



where 

C(5,b) 
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(ii) For any probability measure \x on R, and any 6 > 0, 
F N (d w (jl»,n)-E N [d w (fi%,Li)}>8) 

f NS ( N5 V 2 \ 1 / NR5 2 \ 1 _ 

^^{7^-[^ + ip) h \ 1 + v^)\- (L12) 

Remark 1.6 (i) As in Theorem \1.2\ the dependency ofC(5,b) in 5 and b 
can be made more precise. A key step in the proof of U.ll\) is to choose 
t such that 

E 7V [tr iV (l { |x A |> T} )] < 5/12b, 

and then C(S,b) is determined by r. Minimizing, in t, the right hand 
side of K2. 38\) , leads to the following estimate 

f V 2 /RfMr-Gifrj) 
E Ar [tr iV (l {|XA| > T} )] < exp | - —el ^ y2 L 

where £(x) = (1 + x) ln(l + x) — x. For x > 1, 2£(x) > xlnx. Hence 
one can choose t to be the solution, in x, of the equation 

x xR , 126 
It then follows that C(S, b) can be taken to be 



Outside of the finite exponential moment assumption, an interesting class 
of random matrices with infinitely divisible entries are the ones with stable 
entries, which we now analyze. 

Recall that X in M. d is a-stable, (0 < a < 2), if its Levy measure v is 
given, for any Borel set B G B(M. d ), by 

v{B) = / a(dO / IbK)^, (1-13) 
Js*- 1 Jo r 

where a, the spherical component of the Levy measure, is a finite positive 
measure on S '™ 1 , the unit sphere of lR d . Since the expected value of the 
spectral measure of a matrix with a-stable entries might not exist, we look 
at the deviation from a median. Here is a sample result. 
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Theorem 1.7 Let X = (^i,^j,^j t j)i<i<j<N be an a- stable, < a < 2, 
random vector in M, n2 with Levy measure v given by U.13\) . 

(i) Let f G Lip(l), and let m(trjsr(f(KA})) be any median o/£rjv(/(XA))- 
Then, 

/ oN 2 -l\ 

P JV (tr J v(/(X A ))-m(tr JV (/(X A )))>^)<C( a )(v / 2a)° ^ > , (1.14) 



1 1/" 

a + ea)/a(2 — a). 



whenever SN > y/2* 



2a(^ 2 - 1 )C(a) , and where C{a) = 4 a (2 



(zzj Let \ ma x(X-A) be the largest eigenvalue o/Xa, and let m(X ma x(^A)) 
be any median o/A m „(XA), then 

F N (X max (X A )-m(X max (X A )) >5)< C(a)(V2a) a ^^J , (1.15) 



l/as 

a + ea)/a(2 — «). 



whenever 5vN > \[2< 



2a(S N2 ~ 1 )C(a) , and where C{a) = 4 a (2 - 



Remark 1.8 Lei M be a Wigner matrix whose entries Mi$, l<i<N , Mf 
l<i<j<N, and Mjj, l<i<j<N, are iid random variables, such that the 
distribution of (Mi^l belongs to the domain of attraction of an a-stable dis- 
tribution, i.e., for any 5 > 0, 

P(|Mi 4 | > S) = 

for some slowly varying positive function L such that lim L(tS) / L(5) = 1, 

S— >oo 



for all t > 0. Soshnikov J28y showed that, for any 5 > 0, 

lim ¥ N (X max (b N 1 M) >5) = 1- exp(-r Q ), 

N^oo 

where is a normalizing factor such that lim N 2 L(b at) / 'b% = 2 and where 

N^oo 

X ma x{bjf M) is the largest eigenvalue ofbjj M. In fact lim N~~ e /b^ = and 
lim b]y/Na +e = 0, for any e > 0. As stated in fTEj . when the random vector 
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X is in the domain of attraction of an a-stable distribution, concentration 
inequalities similar to fli.i^l ) or U.15}) can be obtained for general Lipschitz 
function. In particular, if the Levy measure of X is given by 

"(B)= [ a(dO r°°lB(rO^^, (1-16) 

for some slowly varying function L on [0, +oo) ; and if we still choose the 
normalizing factor bjy such that liniAr^oo a(S N ^ 1 )L(bN)/b% is constant, then, 



P JV (A moa; (6^ 1 M) - m(A ma!E (&^M)) > 5) 

^ C{a)a{S N2 -^ L {^ 2/ 

b a N 5* ' 1 J 

whenever 

(5b N ) a > 2 l+a ' 2 C(a)a(S N2 - l )L{b N 5/V2). 

Now, recall that for an N 2 dimensional vector with iid entries, a(S ! 
iV 2 (<j(l) + <r(— where <r(l) is short for cr(l,0, . . . , 0) and similarly for 
<r(— 1). Thus, for fixed N, our result gives the correct order of the upper 
bound for large values of 5, since for 5 > 1, 



e-1 

ed a 



< 1 - e- 5 ~ a < 4- 



Moreover, in the stable case, L(5) becomes constant, and bjq = N 2 l a . Since 
^max{N~ 2 l a Wi) is a Lipschitz function of the entries of the matrix M with 
Lipschitz constant at most \/2N~ 2 / a , for any median m(A ma;c (iV _2 / a M)) of 
\ma X (N- 2 / a M), we have, 

F N (X ma ^N-iM)-m(X max (N-iM))>5)<C(a) ^±^1} 1 (1.18) 



whenever 5> [2C{a)(a{\) +a(-l))] 1/a . Furthermore, using Theorem 1 
[13^ . it is not difficult to see that m(A max (iV~ 2//a M)) can be upper and lower 
bounded independently of N. Finally, an argument as in Remark \1.15\ below 
will give a lower bound on \ max {N~ 2 l a Wi) of the same order as U.18\) . 



in 
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The following proposition will give an estimate on any median of a Lips- 
chitz function of X, where X is a stable vector. It is the version of Proposition 
11.11 for a-stable vectors. 

Proposition 1.9 Let X = {^i,^j,^ij)i<i<j<N be an a-stable, < a < 2, 
random vector in M n2 with Levy measure v given by M.13\) . Let f G Lip(l), 
then 

(i) any median m(f (X)) of f(X) satisfies 

\m(f(X))-f(0)\ 

^■=(^P) ' ( V / i^ +3 <W 1/4) ) +£ ' (" 9 > 

(ii) the mean K N [f(X)} of f(X), if it exists, satisfies 

|E N [/(X)] - /(0)| 

<^)-(^) ' (y=^ + ^(l/4)) + ^, (1.20) 



where k a ui2- a ){x), x> 0, is the solution, in y, of the equation 




(1.21) 



with ei, e2, . . . , e^2 being the canonical basis ofM. . 

Remark 1.10 (%) When the components of X are independent, a direct 
computation shows that, up to a constant, E in both J±(a) and ^(cn) 

is dominated by ( 1 4a — - J , as N — > oo. 
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In complete similarity to the finite exponential moments case, we can 
obtain concentration results for the spectral measure of matrices with a- 
stable entries. 

Theorem 1.11 Let X = (ufi, uifp ^u)i<i<j<N be an a-stable, < a < 2, 
random vector in M. N with Levy measure v given by U.lty) . 

(i) Then, 



P JV ( sup \tr N (f(X A ))-E N itr N (f(X A ))}\ >5 

yf£Lip b (l) 



n a rr( Q^- 1 ) 

<C{8M J a§a j Al, (1.22) 



where 



., , , , , ,' v / 2a\ 1+ VJi(«) + l \ 1+a , 
C(5,b,a)=[C 1 (a)[^=J f— ^ + bj +C 2 (a, 



with Ci(a) and C^a) constants depending only on a, and with J\{a) 
as in Proposition \1.9l 

(ii) For any probability measure fi, 

I oN' 2 -l\ 

F N (d w (^^)-m(d w (^^))>5)<C( a )(V2ar ^ \ (1.23) 



whenever 5N > \[2a 
a + ea)/a(2 — a). 



2a(S N2 ~ 1 )C(a)] 1/a and where C(a) = A a (2 



It is also possible to obtain concentration results for smaller values of 5. 
The lower and intermediate range for the stable deviation obtained in [1] 
provide the appropriate tools to achieve the following result. We refer to |1] 
for complete arguments, and only provide below a sample result. 

Theorem 1.12 Let X = ^/l-, j)i<t<i<N be an a-stable, 1 < a < 2, 
random vector in M n2 with Levy measure v given by U.13\) . For any e > 
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0, there exists i](e), and constants D\ = D 1 (a,a, N,cr(S N2 x )) and D 2 = 
D 2 (a, a, N, a(S N2 ~ 1 )) , such that for all < 8 < 77(e), 

F N ( sup \tr N (f(X A ))-E N [tr N (f(X A ))]\ > 5) 

feLi Pb {i) 

<(l + e)^exp(-D 2 8 2 -^-). (1.24) 



Remark 1.13 (y In ( 1.14 ), /ll-15\) or $1.25\) . the constant C\a) is not 



of the right order as a — ► 2. It is, however, a simple matter to adapt 
Theorem 2 of fWj to obtain, at the price of worsening the range of 
validity of the concentration inequalities, the right order in the constants 
as a — > 2. 

(ii) Let us now provide some estimation of Di and D 2 , which are needed for 
comparison with the GUE results of J^j (see (Hi) below). Let C(a) = 

2 a (ea + 2 - a)/(2(2 - a)), K(a) = max |2 Q /(a - 1), C(a)}, L(a) = 
((a - l)/a) a/{a ~ l) {2 - a)/10 and let 



+ ^(l2C(a)a(^ 2 - 1 ))"6i (1.25) 
v4s shown in the proof of the theorem, Di = 24D* , while 
D 2 = L ^ ( N ) ^ 1 

( \ 1/a 

Thus, as N — > +oo ; D\ is of order N~ 1 ^ 2 I a(S N _1 ) J , while D 2 is 



of order N 3a ^ 2a ~ 2 ^ a^ 2 " 1 ) 
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(Hi) Guionnet and Zeitouni f3j[, obtained concentration results for the spec- 
tral measure of matrices with independent entries, which are either com- 
pactly supported or satisfy a logarithmic Sobolev inequality. In partic- 
ular for the elements of the GUE, their upper bound of concentration 
for the spectral measure is 

i^^r^M^)' (L26) 

where C\ and C2 are universal constants. In Theorem \1.12\ the order, 
in b, of Di is at most b a+l l a , while that of D2 is at least 5~( a + 1 )/(°~ 1 ) . 
This order is thus consistent with the one in U.26]) . as a is close to 2. 
Taking into account part (ii) above, the order of the constants in U.24\ ) 



are correct when a — > 2. Following ^ (see also Remark 4 in P^j)- we 
can recover a suboptimal Gaussian result by considering a particular 
stable random vector and letting a — > 2. Toward this end, let 
X^ a > be the stable random vector whose Levy measure has for spherical 
component a, the uniform measure with total mass o-(S N _1 ) = N 2 (2 — 
a). As a converges to 2, X^ converges in distribution to a standard 
normal random vector. Also, as a — > 2, the range in 5 in Theorem 
\1.12\ becomes (0, +00) while the constants in the concentration bound 
do converge. Thus, the right hand side of U.24\ ) becomes 



which is of the same order, in 5, as U.26]) . However our order in N is 
suboptimal. 

(iv) In the proof of Theorem \1.12\ the desired estimate in 1(2.50) is achieved 
through a truncation of order 5~ l l a , which, when a — > 2, is of the same 
order as the one used in obtaining U.26]) . However, for the GUE result, 
using Gaussian concentration, a truncation of order a/1ii(126 / 5) gives 
a slightly better bound, namely, 



■ exp 



6 " r 1 8ca 2 ln±f 



where C\ and C2 are absolute constants (different from those of \1.26) ). 
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Wishart matrices are of interest in many contexts, in particular as the 
sample covariance matrix in statistics. Recall that M = Y*Y is called a 
complex Wishart matrix if Y is a K x iV matrix, K > N, with entries Yjj = 



variables with finite variance a 2 , the empirical distribution of the eigenvalues 
of Y*Y/iV converges as K — > oo, N — > oo, and K/N — ► 7 G (0, +00) to the 
Mar cenko- Past ur law ([3], |21j ) with density 



where C\ = cr 2 (l — 7~ 1 ^ 2 ) 2 and c 2 = o- 2 (l + 7~ 1 / 2 ) 2 . When the entries of Y 
are iid Gaussian, Johansson [Hj and Johnstone [15] showed, in the complex 
and real case respectively, that the properly normalized largest eigenvalue 
converges in distribution to the Tracy- Widom law ([30], [31]). Soshnikov 
[27] extended the result of Johnstone to Wishart matrix with Non-Gaussian 
entries under the condition that K — N = 0(N 1 ^ 3 ) and that the moments 
of the entries do not grow too fast. Soshnikov and Fyodorov [29] recently 
studied the distribution of the largest eigenvalue of the Wishart matrix Y* Y, 
when the entries of Y are iid Cauchy random variables. We are interested 
here in concentration for the linear statistics of the spectral measure and 
for the largest eigenvalue of the Wishart matrix Y*Y, where the entries of 
Y form an infinitely divisible and, in particular, a stable one. We restrict 
our work to the complex framework, the real framework being essentially the 
same. 

It is not difficult to see that if Y has iid Gaussian entries, Y*Y has 
infinitely divisible entries, each with a Levy measure without a known explicit 
form. However the dependence structure among the entries of Y* Y prevents 
the vector of entries to be, itself, infinitely divisible (this is a well known fact 
originating with Levy, see |25j). The methodology we previously used cannot 
be directly applied to deal with functions of eigenvalues of Y*Y. However, 
concentration results can be obtained when we consider the following facts, 
due to Guionnet and Zeitouni [9] and already used for that purpose in their 
paper. 




p 7 {x) 



1 



a/ (c 2 - x)(x - Ci), Ci < X < c 2 , 



2lTX r )<J 2 
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Let 



and 



hJ 



then Xa = 




for 1 < % < K, 1 < j < K 

for N+l<i<K + N,K + l<j<K + N 

forl<i<K,K+l<j<K + N 

for N + 1 < i < K + N, 1 < j < K, 



iorl<i<K,l<j<K 
iorN+l<i<K + N,K + l<j<K + N 
iorl<i<K,K+l<j<K + N 
for N + 1 < i < K + N, 1 < j < K, 



(1.27) 



(1.28) 



G -M(^ +A r) x (^+AT)(C), and 



xi = 



Y*Y 







YY* 



Moreover, since the spectrum of Y*Y differs from that of YY* only by the 
multiplicity of the zero eigenvalue, for any function /, one has 



and 



tr(/(Xi)) = 2tr(f(Y*Y)) + (K — N)f(0), 



X max (M 1/2 ) = max |A,(X A )|, 



where M 1 / 2 is the unique positive semi-definite square root of M = Y*Y. 

Next let ¥ K ' N be the joint law of (Yg, Yj^i^K^jKN on R 2KN , and let 
E K ' N be the corresponding expectation. We present below, in the infinitely 
divisible concentration result for the largest eigenvalue X max (Wl), of 

the Wishart matrices M = Y*Y. The concentration for the linear statistic 
£ r Jv(/(M)) could also be obtained using the above observations. 
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Corollary 1.14 Let M = Y*Y , with Y id = Yg + V^IY 1 ^. 

(i) Let X = (Yfj, j)i<i<N,i<j<K be a random vector with joint law 
f K > N ~ ID(/3,0,v) such that E K ' N [e t ^] < +00, for some t > 0. 
Let T = sup{t > : E K < N [e tm ] < +00} and let h' 1 be the inverse of 

h(s)= I |M|(e slHI - l)v(du), 0<s<T. 

Then, 

V K > N (A_(M^) - E K ' N [X max (M^)] >5)< e-tf'^M* (1.29) 
for allO < 5 < h{T~). 

(ii) Let X = (Yij, ^lj)i<i<K,i<j<N be an a-stable random vector with Levy 
measure v given by v{B) = J s2K n-i cr{d£) / + °° ls(r£)<ir /r 1+a . Then, 

P^A^M 1 / 2 ) - m(X max (M 1 / 2 )) >5)< C(a)(V2) a ^-^ L , 

whenever 5 > V2a [2a(S 2KN ~ 1 )C( y a)] 1/a and where C(a) = 4 a (2-a + 
ea)/a(2 — a). 

Remark 1.15 (i) As already mentioned, Soshnikov and Fyodorov (129^) 
studied the asymptotic for the largest singular value of the K x N ran- 
dom matrix Y , which is the largest eigenvalue of the Wishart matrix 
Y*Y, when the entries of Y are iid Cauchy random variables. They 
argue that although the typical eigenvalues of Y*Y is of the order KN , 
the correct order of the largest eigenvalue of such a matrix is K 2 N 2 . 
Our result implies that the largest eigenvalue A max (M) of the Wishart 
matrix M = Y*Y, when the entries of Y form an a-stable random 
vector, is of order at most a(S 2KN ~ 1 ) 2 ^ a . We also have a lower bound 
result which is described next. In particular, if the entries of the matrix 
Y are iid a-stable random variables, the largest eigenvalue of Y*Y is 
of order K 2 l a N 2 ' a . 

(ii) Let X ~ ID(P, 0, v) in ~R d , then (see Lemma 5.4 in f^) for any x > 0, 
and any norm \\ ■ ||jv on M. d , 

F(\\Xy>x) > ^(l-exp{-v({ueR d : \\u\\ M > 2x})^. 
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But, A max (M 1//2 ) is a norm of the vector X = (Y^-,Yf •), which we 
denote by \\X\\\, if X is a stable vector in M. 2KN . 



^(X max (MV*)- m (X max (M^)>6 
= P^A^M 1 / 2 ) > 5 + m(A ma:E (M 1 / 2 )) 
> ^(l-exp{ - u({X max (M 1 ^ 2 ) > 2(5 + m(A ma:c (M 1 / 2 )))})} 
I (i _ eX p { - u({ ||X|| A > 2(5 + miXma^M 1 / 2 ))) }) } 



> 

~ 4 

= H 1_eXP l a ( 5 + m (A m I(MV 2 )))^J • ^ L3 °) 

where S 2 ^ -1 zs i/ie um'£ sphere relative to the norm || • ||a and where 
a is the spherical part of the Levy measure corresponding to this norm. 
Moreover, if the components of X are independent, in which case the 
Levy measure is supported on the axes ofM. 2KN , a^S 2 ^ -1 ) is of order 
KN , and so the largest eigenvalue o/M 1 / 2 is of order K l l a N l l a . 

(Hi) For any function f such that g(x) = f(x 2 ) is Lipschitz with Lips- 
chitz constant \\g\\ L ip ■= \\\f\\\c, tr(g(X A )) = tr(f(X. 2 A )) is a Lip- 
schitz function of the entries of Y with Lipschitz constant at most 
\^2\\\f\\\cVK + N. Hence, under the assumptions of part (i) of Corol- 
lary ^. 14 , 



F K ' N (tr N (f(M)) — E K,N [tr N (f(M))] > 8* 1 A 



N 

y/2{K+N)s/\\\f\\\c 

<exp<{ - / hr l {s)ds\, (1.31) 







for all < 6 < \\\f\\\ c h (T~) fy/2(K + N). 



(iv) Under the assumptions of part (ii) of Corollary 1.14, for any function 



f such that g(x) = f(x) is Lipschitz with \\g\\up — \ \\f\\\c, o,ny median 
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m(tr N (f(M))) oftr N (f(M)), 

.K + N 



P K,iV ( tr N (f(M)) - m(tr N (f(M))) > 5- 



N 



la fT (c2KN-l\ 

<C(a) . "unl^ ^ J (132) 

w/ienewer 5 > 1 1 1/| 1 1 £ ^(S 2 *"- 1 )^)] 1/a / ^2(kTN), and where 
(7(a) = 4 Q (2 - a + ea)/a(2 - a). 

Remark 1.16 TTie methodology used to obtain the results of the present pa- 
per, in the absence of the finite exponential moments, can be applied to any 
matrices whose entries on and above the main diagonal form such an in- 
finitely divisible vector X . However, to obtain explicit estimates, we do need 
specific bounds on V 2 {r) and v(r), which are not always available when fur- 
ther knowledge on the Levy measure of X is lacking. 



2 Proofs: 

We start with a proposition, which is a direct consequence of the concentra- 
tion inequalities obtained in [11] for general Lipschitz function of infinitely 
divisible random vectors with finite exponential moment. 

Proposition 2.1 Let X = ^/j)i<i<j<iv be a random vector with 

joint law F N ~ ID((3, 0, v) such that E N [e* l|x|1 ] < +oo ; for some t > and 
let T = sup{t > : E^e*^] < +00} . Let h' 1 be the inverse of 

h(s) = [ \\u\\ (e slHI - l)v{du), < s < T. 
Jr n2 

(i) For any Lipschitz function f, 

-J^ amLiv h- 1 {s)ds 

for all < 8 < V2a\\f\\ Up h (T')/N. 
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(ii) Let Xmax(^-A) be the largest eigenvalue of the matrix X A - Then, 

F N (X max (X A ) - E N [X max (X A )] > 5) < exp | - br\s)dA, 
for allO < 5 < V^ah (T~) /^fN. 



Proof of Theorem 11.2b 

For part (i), following the proof of Theorem 1.3 of PJ, without loss of 
generality, by shift invariance, assume that min{x : x G JC} = 0. Next, for 
any v > 0, let 

{0 if x < 
x if < x < v (2.33) 
v if x > v . 

Clearly g v G Lip{l) with ||<7t,||oo — u - Next for any function / G Lipic(l), any 
A > 0, define recursively f A (x) = for x < 0, and for (j — 1)A < x < jA, 
j = l,...Jfl,let 

r— i 

/a(s)=X>2\ 
3=1 

where g^' := (21{/(,-A)>/ A ((,-i)A)} - 1)0a(s - (j - 1)A). Then \ f - f A \ < A 
and the 1-Lipschitz function /a is the sum of at most \IC\/ A functions g^ G 
Lip(l), regardless of the function /. Now, for 6 > 2A, 



( 



¥ N sup \tr N (f(X A )) - E N [tr N (f(X A ))}\ > 6 
\feLi PK (i) 



<F N ( sup {\tr N (f A (X A ))-E N (tr N (f A (X A )))\ + \tr N (f(X A )) 
\feU PK (i) I 

-tr N {f A (X A ))\ + \E N [tr N {f(X A ))} -E N [tr N (f A (X A ))}\ } > 5 



< P JV swp\tr N (f A (X A )) - E N (tr N (f A (X A )))\ >5-2A 

V /a 
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< M sup ( \tr N (gf(X A )) - E N [tr N (g^(X A ))]\ > 



^2.34) 



whenever < 5 < y Sv^ol/CI/i (T ) /iV, and where the last inequality follows 
from part (i) of the previous proposition by taking also A = 5/4. 

In order to prove part (ii), for any / G Lip b (l), i.e, such that ||/|| Lip < 1, 
, < b, and any r > 0, let f T be given via: 



/rO) = < 



'/(^) if M < r 

/(r) - sign(/(r))(x -r) if r < x < r + |/(r)| 

/(-r) + sign(/(-r))(x + r) if -r - \f(-r)\ < x < -r 

otherwise. 



(2.35) 



sup 



< sup 



< sup 



Clearly f T e Lip(l) and supp(f T ) C [— r — |/(— r)|,r + |/(t)|]. Moreover 
tr iY (/(X A ))-E JV (tr JV (/(X A ))) 

tr JV (/ r (X A ))-E iV (tr i v(/ r (X A ))) i 

+ sup tr iY (/(X A )-/ r (X A ))-E Ar [tr JV (/(X A )-/ r (X A ))] 

feLi Pb (i) 

tr N (f T (X A )) - E N (tr N (f T (X A ))) 
+ 2tr N (g b (\X A \ - r)) + 2E A > JV (^(|X A | - r))], (2.36) 
with g fo given as in (I2.33p . Now, 



F N ( sup \tr N (f(X A ))-E N (tr N (f(X A )))\>5 

\f£Lip b (l) 



<F N ( sup \tr N (f T (X A ))-E N (tr N (f T (X A )))\> d - 
VeLiPb(i) d/ 

+ P w (2tr N (g b (\X A \ - r)) + 2E N [tr N (g b (\X A \ 



r))}> 



25^ 
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<W sup \tr N (f T (X A ))-E N (tr N (f T (X A )))\> 6 - 



+P^(^(|X^-r))-E> JV (^(|X A |-r))]>^-2E> w (^(|X^-r))] 



<W sup |^(/ T (X A ))-E Ar (tr iV (/ T (X A )))|>- 
v /eLi Pi) (i) «5 

+P 7V (tr iV (^(|X A |-r))-E JV [tr iV (^(|X A |-r))]>--26E Ar [tr Ar (l { |x A |> r 



^2.37) 



Let us first bound the second probability in (12.371) . Recall that the spectral 
radius p(X A ) = max |Aj| is a Lipschitz function of X with Lipschitz constant 

at most a^j2/N. Hence, for any < t < T, and 7 > such that i>(p 7 ) < 1/4, 

1 N 



1=1 



<P JV (p(X A )>r) 

< exp U(t) - (^-t - G 3 ( 7 )V1 (2.38) 

where we have used Proposition 11.11 in the next to last inequality and where 
the last inequality follows from Theorem 1 in [TT] (p. 1233) with 

H(t)= [ h{s)ds= [ (e t|MI -t\\u\\ - l)v{du). 
Jo Jr n ' 2 

We want to choose r, such that E N [tr N (l{\x_ A \>r})] < 5/126. This can be 
achieved if 

VN „ , , In + H{t) , 

r-G 2 ( 7 )> s — — . (2.39) 



V2a t 

Since 

''1 



d f\n If + H(t) \ _ th(t) - In ^ - H(t) 



dt\ t J t 2 
24 



and 

d 2 flnif + H(t)\ t 3 H"{t) -2t{th{t) -In ^ -#(*)) 



dt 2 \ t J t A 

it is clear that the right hand side of (I2.39f) is minimized when t = to, where 
to is the solution of 

th(t)-H(t)-\n^ = 0, 
o 

and the minimum is then h(to). 
Thus, if 

r = Co(6, b) := ^ ( G 2 ( 7 ) + h(to)) , (2.40) 



/iV 

then 



E JV [tr^(l { | XA |>.})] < ^, 



and so, 



P W (^(^(|X^-r))-E> iV (( 7fe (|X^-r))]>^-26E JV [tr w (l { |x A |> T} 
<P iY ^(^(|X A |-r))-E JV [tr JV (^(|X A |-r))]>^ 

< exp | - jf 6v ^ a /T 1 (s)ds | , (2.41) 

for all < 5 < 6\/2ah (T~) /N, where Proposition 12.11 is used in the last 
inequality. 

For t chosen as in (12.401) . let /C = [— r — b, r + b], then for any / e Lipi>(l), 
f T G Lipic(l). By part (i), the first term in ( 12. 37ft is such that 



F N ( sup |tr w (/ T (X A ))-E Ar (tr^(/ T (X A )))| > i 
y feLi Pb (i) ■> 



< P JV sup |tr^(/ T (X A )) - E N itr N (f T (X A ))]\ > 
^48(C ( < S ; 6) + &) exp f _ r^Mc oi s, b)+ » h _ 1{s)ds . 



5 

for all 0<5 2 < 144v^a(C (5, 6) +6)/i(T-) /iV. 
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Hence, returning to (12.371) . using 02.411) and f!2.42p and for 



5 < min \ 6V2ah (T~) /N,JluV2a(C (5, b) + b)h(T~)/N 



we have 

F N ( sup \tr N (f(X A ))-E N (tr N (f(X A )))\>5 
\feiAp b (i) 



<| 2+ i.) gim±j) exp 



12 



144 V2a(C (5,b)+b) 



h~ l {s)ds}, (2.43) 



since only the case 5 < 2b presents some interest (otherwise the probability 
in the statement of the theorem is zero). Part (ii) is then proved. 

□ 



Proof of Proposition II. 4t 

As a function of x G M. N , d^lfi 1 ^, fi)(x) is Lipschitz with Lipschitz con- 
stant at most \^2a/N. Indeed, for x, y G M. N , 



dw(fi>A>t i )( x )= SU P 
feLi Pb (i) 

< sup 

feLi Pb (i) 



tr N (f(X A )(x)) - [ fd/i 
Jw 

tr N (f(X A )(x))-tr N (f(X A )(y)) 



+ sup 



tr N (f(X A )(y))- / fdfi 



Theorem 11.41 then follows from Theorem 1 in [TT] . 



(2.44) 
□ 



Proof of Corollary II. 5t 

For Levy measures with bounded support, ~E N [e*" x "j < +oo, for all t > 0, 
and moreover 

hit) < V 21 



e tR -l 



R 
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Hence ^ 

(t) = J h{s)ds < ^ (s ti? -1-tR), 

and 

exp { - jTVwd.} < exp {I - (| + I!) in (l + }. 
Thus, one can take 

C{6 > h) = C ( W ( Gs (7) + T {et ° R " 1} ) + b ) ' 

where to is the solution, in t, of 

Applying Theorem 11.21 (ii) yields the result. 

□ 



In order to prove Theorem II. Ill we first need the following lemma, whose 
proof is essentially as the proof of Theorem 1 in [12] . 

Lemma 2.2 Let X = (co^, u>^, ^lj)i<i<j<N be an a-stable vector, < a < 
2, with Levy measure v given by M.13\) . For any x ,Xi > 0, let g xo>xl (x) — 
g X i{ x ~ x o)> where g xi (x)is defined as in 112.33]) . Then, 



P 



N 



tr N (g X0 , Xl (X A ))-E N [tr N (g X(hXl (X A ))} > S < C(a 



a a a(S 



AT2_1N 



N a S a 



whenever 5 1+a > (2V2~a) +a a{S N1 )xi/aN 1+a and where C {a) = 2 5a / 2 (2ea+ 
2-a)/a(2-a). 



Proof of Theorem 11.111 

For part (i), first consider / G Lip;c(l). Using the same approximation 
as in Theorem 11.21 an y function / G Lip K {l) can be approximated by /a, 
which is the sum of at most |/C|/A functions g^ G Lip(l), regardless of the 
function /. Now, and as before, for 5 > 2A, 
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P"[ sup \tr N (f(X A ))-E N (tr N (f(X A )))\>5 



< ^ sup P J 



fr^(XA))-E w [tr^(X A ))] 



> 



|/C| 



^4|^| 8»a»C 2 ( a )g(^ 2 - 1 )|/Cr ,„ ... 

whenever 

* 2 >^f^)- (2 . 46) 



8|/C| V 4a 

and where the last inequality follows from Lemma [2.21 taking also A = 5/4. 

For any / £ Lipt,(l), and any r > 0, let f T be given as in (12. 35ft . Then, 
f T £ Lipic(l), where /C = [— r — b, r + b], and moreover, 



f sup \tr N (f(X A ))-E N (tr N (f(X A )))\ > s] 
\feLi Pb (i) J 



<P^(tr^(^, b (|X^))-E JV [tr iV (( ?Ti6 (|X^))]>--2ffi Ar [tr J v(l { |x A |> T 

/7'v(./r(X A )) - E-M/r.v(/ r (X A )))| > ^ 

K freLi PK (l) 



+ F N ( sup |^(/ r (X A ))-E Ar (tr w (/ T (X A )))|>-). (2.47) 



The spectral radius p(X A ) is a Lipschitz function of X with Lipschitz con- 
stant at most \/2a/vN. Then by Theorem 1 in [T2] . 



i=l 

<P 7V (p(X A )> 



< P^ ( p(X A ) - m(p(X A )) > r - ^£ Ji(a) 



AT 



< C 1 (a)2°/ 2 aV(^- 1 ) 
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whenever 

A^)-,?*)^, (2 . 49) 



Hsi J ~ N a l 2 
and where C±(a) = 4 a (2 — a + ea)/a(2 — a). Now, if r is chosen such that 

C l {a)2 a ' 2 a a a{S N2 - 1 ) 5 



: J M) > TTTZn , (2.50) 



that is, if 



it then follows that 

Since g T ,b(\^-A) * s the sum of two functions of the type studied in Lemma [2T21 
with x\ = b, we have, 



P JV (tr iV (^(|X A |))-E iV [tr w (^ 6 (|X A |))]>--26E 7V [^(l { |x A |> T} 
< 2¥ N (tr N (g T)b (K A )) - E N [tr N (g T;b (X A ))} > — 

£2C2(a) H!^p), ( , 51) 



S > + * > f^'^l^M, (2 .52) 



whenever 

v iv y a 

and where C^a) = 2 5a ^ 2 (2ea + 2 — a)/a{2 — a). The respective range 
and (I2.52p suggest that one can choose, for example 

V2a \[2a 
t = —=Ji{a) H =6. 

Vn Vn 
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Then, there exists 6 (a, a, N, v) such that for 5 > 8(a, a, N, u), 



P JV sup |tr w (/(X A ))-E JV [tr JV (/(X A ))]| >5 

v /eLi Pi) (i) 

<F N ( sup \tr N (f T (X A ))-E N [tr N (f T (K A ))]\ > 

V t l; pk (i) 



+ P w (tr^(^, fe (|X A |))-E w [tr^(^(|X A |))]>--26E Ar [tr iV (l { |x A |> r} ] 

l+Q 



< C 3 (a)aV(^ 2 - 1 )(^Ji(«) + 6+^ / 



where C 3 {a) = 2 4+2a 12 a C 2 (a), C 4 (a) = 2(12 a )C 2 (a) and 5(a, a, N, v) is such 
that fl2^6|) and f[2~52D hold. 

Part (ii) is a direct consequence of Theorem 1 of [T2], since dwif^A^) 
Lip(\r2a/N) as shown in the proof of Proposition II .41 □ 

Proof of Theorem 11.121 

For any / e Lip(l), Theorem 1 in [12] gives a concentration inequality for 
f{X), when it deviates from one of its medians. For 1 < a < 2, a completely 
similar (even simpler) argument gives the following result, 

P N (f(X) -E N [f(X)]>x) < ^If^ , (2.53) 

whenever x a >K(a)a(S N2 ~ 1 ), where C(a) = 2 a (ea + 2 — a)/(a(2 — a)) and 
AT(a) = max {2 a /(a - 1), C(a)}. 

Next, following the proof of Theorem 11.21 approximate any function / e 
Lipb(l) by / T e Lip[_ T _b )T+ 6](l) defined via (j2.35|) . Hence, 



W sup |tr JV (/(X A ))-E 7V (tr JV (/(X A )))| > <j) 
v /eLi Pb (i) y 

<P^( sup |tr JV (/ r (X A ))-E JV (tr^(/ T (X A )))|>^ 

r 

+P^(tr^(^(|X A |))-E iV [tr^(^ 6 (|X A |))]>--26E Ar [tr J v(l { |x A |> T }]). 



(2.54) 
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For p(Xa) the spectral radius of the matrix Xa, and for any r, such that 
r-E N [p(X A )} > ^K(a)a(S N2 -i)) 1/a , 

E N (tr N (l {lXAl>T} )) < P 7V (p(X A ) - E N \p(X A )] > r - E"[p(X A )]) 

~ (r-E^[p(X A )]) Q ' (2 - 55) 

where we have used, in the last inequality, (I2.53j) and the fact that p(Xa) G 
Lip(^k). For Q > 0, let r = E Af [p(X A )] + With this choice, we 

then have: 



(^] a C{a)a{S N2 - 1 ) 
E-(^(W. } ))< ^ E ^ (XA)]r 

(^)W- 1 ) 

< (2.56) 
- 126' v 7 

provided Q a /5 > V^aK^aiS^-^/VN, and C{a)a(S N2 ' 1 ) / Q a < 

1/(126). Now, taking Q = yfa(mC(a)o(S^- 1 )) 1,a /y/N, and recalling, for 
1 < a < 2, the lower range concentration result for stable vectors (Theorem 
1 and Remark 3 in [3]): For any e > 0, there exists i] (e), such that for all 
< 5 < V2a\\f\\ LipVo (e)/N, 



P"(tM/(X A )) -E N (tr N (f(X A ))) > 6) 

2-a ( a-l\^=- 
10 V a J 

(a(S" 2 -i)y/(«-V\V2a\\f ULlii 



2-a I a-l \ q_i / AT \ "~ 1 

<(l + eW - /° [*> [ * 5 



Q - 1 



(2.57) 

With arguments as in the proof of Theorem 11.21 if 
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there exist constants Di(a, a, N, <t(S n2 1 )) and D 2 (a, a, N, a(S N2 1 )), such 
that the first term in ( 12.541) is bounded above by 



Di(a,a : N,a(S N2 ~ 1 )) 



o a 



exp(-D 2 (a, a, N, aiS^'^S^) . (2.58) 



Indeed, with the choice of r above and D* as in fll.25p . 2(r + b) < D*/5 1 ^ a . 
Moreover, as in obtaining (I2.34p . Di can be chosen to be 24/}*, while D 2 can 
be chosen to be 

a 

2— a I a—1 \ a-i / at 
10 \ a ) ly 



^(S* 3 - 1 ))^^ ' ( 72D *) 



We remind the reader that, as already mentioned, J^ac) can be replaced 
by E [||X||]. According to the result of Marcus and Rosihski [20] and the 
estimate in [13], if E N [X] = 0, then 



a{s , -x } < E ,v [||x||] < _ 



4(2 -«) x / a ^ 8((2-a)(a-l)) 1/a 

Finally, note that, as in the proof of Theorem 11.21 (ii), the second term in 
( I2.54p is dominated by the first term. The theorem is then proved, with the 
constant Di(a, N,a(S N _1 )) magnified by 2. 

□ 

Proof of Corollary I1.14t 

As a function of (Y^-, j)i<i<K,i<j<N, with the choice of A made in 
f ll.27p . \ r nax(X.A) G Lip(y/2). Hence part(i) is a direct application of Theo- 
rem 1 in pTJ, while part(ii) can be obtained by applying Theorem 11.71 □ 
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